Effects of Short-Range Correlations in (e, e'p) 
reactions and nuclear overlap functions 

Q ■ M. Mazziotta, J. E. Amaro and F. Arias de Saavedra 

o" 

(N 

Departamento de Fisica Moderna, Universidad de Granada, Granada 18071, Spain 
^' 

\D 

> 

oo 

Abstract 

o 

^^ , A study of the effects of short-range correlations over the (e, e'p) reaction for low missing en- 

ergy in closed shell nuclei is presented. We use correlated, quasi-hole overlap functions extracted 
from the asymptotic behavior of the one-body density matrix, containing central correlations of 
Jastrow type, up to first-order in a cluster expansion, and computed in the very high asymptotic 
region, up to 100 fm. The method to extract the overlap functions is checked in a simple shell 
model, where the exact results are known. We find that the single-particle wave functions of 
Slj . the valence shells are shifted to the right due to the short-range repulsion by the nuclear core. 

The corresponding spectroscopic factors are reduced only a few percent with respect to the shell 
model. However, the (e, e'p) response functions and cross sections are enhanced in the region of 
the maximum of the missing momentum distribution due to short-range correlations. 
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1 Introduction 

In the last years, quasi-free (e, e'p) reactions have proved to be an ideal tool to study 
the spectral function of nuclei |1|, |^ . Apart from the ambiguities coming from final-state 
interactions and off-shell effects in the electromagnetic current, it is possible to extract 
such valuable information about single-particle properties as momentum distributions and 
spectroscopic factors of the different nucleon shells near the Fermi level [^. 

The observed values of the spectroscopic factors for closed shell nuclei are around 
S ~ 0.6-0.7 @|-@] — yet theoretical studies appear to indicate that these values increase 
due to relativistic effects . The small values of the spectroscopic factors are attributed 
to the departure from the extreme single-particle model, where the nucleon is ejected from 
a well defined orbit within the target nucleus. 

Recent Variational Monte Carlo calculations including short-range correlations in the 
nuclear wave function report spectroscopic factors differing from unity only a few percent 



H^Ell- I^ addition, center of mass correlations can enhance the spectroscopic factor by 
a factor around 7% for the valence shell of ^^O |T5|, |I6|. Although it has been shown that 



these values could be lowered by the inclusion of low-energy configuration mixing in the 
wave functions |]T3, the experimental spectroscopic factors, and hence the cross section 



of (e, e'p) reactions for low missing energy, are not satisfactorily explained by present 
theoretical models. 

It has been shown recently ||l3 that the overlap functions between the nuclear ground 



state with A nucleons and a residual state with A — 1 nucleons — which are a main ingredi- 
ent in nucleon knock-out calculations, including implicitly the spectroscopic factor — can 
be obtained from the asymptotic behavior of the one-body density matrix (OBDM). This 
was confirmed in an exactly solvable many-body system in one dimension with a zero- 
range interaction [0. The theorem was applied in |]T^| to obtain the overlap functions 



and spectroscopic factors from a model OBDM including Jastrow correlations, and later 
[^ , |2T[| applied to the calculation of cross sections for several reactions, such as {p,d), 
{e,e'p) and (7,p). However, some problems with the numerical restoration procedure 
performed in |T^ , which is based on an exponential fit of the OBDM, were pointed out in 
p!5[] , namely: i) the use of harmonic-oscillator single-particle wave functions, which do not 
have the correct exponential behavior in the asymptotic region, to construct the OBDM, 
and ii) the nature of the exponential fit allowed to obtain in |TP| overlap functions for the 



unoccupied states, while in |T^ it was demonstrated that such extraction is not possible 
starting from a CBF-tj^e wave function. 



In reference ||21| the same restoration procedure was applied to compute overlap func- 
tions using several correlated OBDM taken from the literature. In particular, the one 
obtained in p2| by a truncated cluster expansion of the radial multipoles pi (r, r') of the 



OBDM. Using the above density, spectroscopic factors corresponding to the pi/2 and ps/2 
shells of ^^O were reported [^ to be identical in this model and equal to 0.981. The 



equality of these values is probably due to the fact that the restoration procedure in |gT 



started from a OBDM for / = 1 which (i) was not separated in spin partners j = 1/2, 3/2, 
and that (ii) was computed up to distances of r,r' = 11 fm, which are not large enough 
to separate asymptotically the two contributions pi/2 and p3/2- This fact, together with 



the problems pointed out above, makes possible that some of the effects attributed in ^1 
to short-range correlations, actually be a consequence of the extraction procedure of the 
overlap functions. 

One of the motivations for this work is to clarify this situation, in particular to explore 
the possibility of separating the two spin-orbit partners j = I ± 1/2 starting from a /- 



multipole of the OBDM. Note that, in the correlated model of p2[, the underlying Slater 
determinant was built with Woods-Saxon single-particle wave functions including spin- 
orbit splitting, so the model can deal with different overlap functions, with different 
energies, corresponding to the pi/2 and lp3/2 shells in ^^O or ^3/2 and ^5/2 in ^°Ca. 

Another aim of this paper is to perform a numerical study of the convergence of 
the asymptotic methods used to extract the overlap functions from the OBDM. We first 
carry out such study in the nuclear shell-model, where the exact overlap functions are 
known (they are equal to the single-particle wave functions). This allows us to determine 
optimum upper values of the coordinates in which one should know the OBDM in order 
to obtain convergence in the extraction procedure. We continue evaluating more realistic 
overlap functions starting from the correlated OBDM of |2^, computed up to ~ 100 fm. 



so that we can check the convergence of the results. In particular, we will be able to 
obtain precise values of the spectroscopic factors for quasi-hole states. 

Finally, in this work the overlap functions resulting from the above task will be inserted 
in a model of the (e, e'p) reaction, in order to evaluate the effects of short-range correlations 
over nuclear response functions and cross sections. We use the DWIA model of [^, ^ 
which includes a new expansion of the relativistic electromagnetic current in powers of 
the missing momentum, but it is not expanded in q or uj. Combined with relativistic 
kinematics, this relativized model was shown in ref. [^ to give the same results as 
the relativistic Fermi gas for the electromagnetic inclusive responses in nuclear matter. 



Moreover, the present model was compared in [26| with a fully relativistic DWIA model 
of the reaction for \Q^\ = 0.8 (GeV/c)^, giving a reasonable description of the Atl 
asymmetries recently measured in ^^O ||27|| . 

This paper is organized as follows. In sect. 2 we summarize the DWIA model we 
use for coincidence electron scattering, and its relation with the overlap functions, with 
some details of the multipole analysis of response functions placed in the appendix. In 
sect. 3 we present the different asymptotic procedures to extract the overlap functions 
from the exact OBDM, and the correlated model of the OBDM including short-range 
correlations of the Jastrow type. In sect. 4 we present a study of the reliability and 
convergence of the asymptotic methods in the shell model. The reader is directed to 
sect. 5 for discussion of the results obtained with the correlated model, where the effects 
of short-range correlations over (e, e'p) observables, overlap functions, and spectroscopic 
factors are analyzed, with a brief application of the model to (7,p) reactions. Finally, our 
conclusions are summarized in sect. 6. 
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Figure 1: One-photon exchange diagram for the (e, e'p) reaction. 

2 DWIA model of (e, e'p) reactions 



2.1 Cross section and response functions 



In this section we summarize those aspects of the (e, e'p) reaction that are of relevance 
for this paper. We consider the process shown in the Feynman diagram of figure 1. Here, 
an incident electron with four-momentum Kj^ = {e, ke) interacts with a target nucleus, 
exchanging a virtual photon with four-momentum given by Q'^ = {Ke — K'^)'^ = {uj,q), 
with K'J^ = {e', kg) the scattered electron four-momentum. The outgoing proton with 
four-momentum P'^ = {E', p') is detected in coincidence with the scattered electron. 

The wave function of the spin- zero target in the ground-state is denoted |$o )) with 
(non relativistic) energy Eq . We are interested in the low-missing energy region, where 
the residual nucleus is left in a bound state, |<l>^^~^)) = |Jq,Mq,), with (non relativistic) 
energy EJ^~^\ We neglect recoil and assume parity conservation. 

We work in the laboratory system, with the z-axis pointing into the q-direction and 
the X-axis in the electron scattering plane. In this reference system, the cross section 
for the (e, e'p) reaction, assuming plane waves and the extreme relativistic limit for the 



electron, can be written in the form p9 



dV 



de'dVt'JVtp, 



hA 



(1) 



where Qp/ = {6', 0') are the proton emission angles, h is the electron helicity, S is the 
unpolarized cross section and A is the electron polarization power. These functions are 
given by 



E = Ka 



M 



vlW^ + vtW^ + vtlW^^ cos (p' + vttW^^ cos 20'J (2) 

A = K a mVt L' W^ ^' sin (j)'. (3) 

Here the kinematical factor K is proportional to the momentum p', 

K=J^^ (4) 

(Tm is the Mott cross section, and vk are factors containing the dependence on the electron 
kinematics: 
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(5) 



Note that in this work the vtl and vtu variables have an extra v2 factor with respect 
to the corresponding definition of ref . 1 2^ . 



The five exclusive nuclear response functions W^ are defined by the following linear 
combinations of longitudinal (L) and/or transverse (T) projections, with respect to the 
transfer momentum q, of the relevant tensor describing the hadronic part of the emission 
mechanism 

cos (j)'w^^ = iy°^ + iy^° cos 2(j)'w^^ = wyy - w'' (6) 

sin0'iy^^' = 2(iy°^-iy^°) 

Here the hadronic tensor, W^'^, is related to the transition matrix elements of the nuclear 
electromagnetic current operator J'^(q). It is defined by 

and it represents the maximum information that can be obtained in these kinds of ex- 
periments. Note that the dependence on the azimuthal angle of the emitted proton, 0', 
is given explicitly in Eqs. @. The final hadronic states entering in the definition of the 
hadronic tensor, \p'ms,^^~^^) are in principle the exact scattering states with the cor- 
responding boundary conditions, i.e., they correspond asymptotically to a nucleon with 
momentum p' and third spin component iris, and a daughter nucleus in the state \^)^~^'). 



2.2 DWIA and overlap functions 

In this paper we consider a DWIA model of the current matrix elements between hadronic 
states that enter into the hadronic tensor (|^. This model is based in the impulse ap- 
proximation, in which we assume that the nuclear electromagnetic current is a one-body 
operator. Hence we are neglecting two-body contributions coming mainly from meson- 
exchange currents. The contribution from these two-body currents was analyzed in [pO[] . 
However, in this work we are not including that contribution, since we are interested in 
studying the genuine short-range correlation effects. 

The impulse approximation current operator can be written in momentum space as 

>(q) = J d'k J^(q + k, k)aj_,kak. (8) 

where J^(q+ k, k) is the single- nucleon current for which we use a new non-relativistic 
expansion to first order in k/M, first proposed in refs. |]25|, ^. The time component 



of this current contains charge and spin-orbit contributions, while the transverse current 
is given as the sum of magnetization plus convection pieces. However, each piece of 
the current differs from the traditional non-relativistic one, containing, in addition to the 
nucleon form factors, relativistic correction factors which depend on q and uj. In ref. p^ it 



was shown that, beginning with the usual non-relativistic Fermi gas and using relativistic 
kinematics plus the new currents, the same longitudinal and transverse inclusive response 
functions are obtained essentially as in the relativistic Fermi gas model for arbitrary values 
of g, which can be bigger than 1 GeV. Note that the usual non-relativistic forms of the 
current, that are also expanded in powers of g/M, begin to fail for high values of g > 500 
MeV/c and cannot be applied for q ~ M. 

The second approximation used in DWIA is the description of the ejected proton 
state as a single-particle wave function x-p' (r) , obtained as the solution of the Schrodinger 
equation with a complex optical potential fitted to elastic-scattering data. The final 
hadronic state is then written as 

l<^L^^'^Xp') = /^'A:Xp'(k)aL|$(^-^)) (9) 

where ^ cPkxpi{\^al_ is the field operator creating a nucleon in the state Xp'(k), i-e., the 
wave function of the ejected proton in momentum space. From here we can write the 
matrix element of the current (|^) as 

(<^L^"'^Xp'l>(q)l$i^^) = /rf^A:ci3A:'x;^(k')J^(q + k,k)(<I.(^-i)|ak.ai+,ak|«l>S^)) (10) 
Now we use the anti-commutation rules of the Fermion operators 

"k'^i+k = '^(k' - q - k) - a]j+k«k'- (11) 

The contribution of the second term in this equation, aq_^j^ak', can be neglected if 



1. the momentum of the ejected proton is much larger than the Fermi momentum of 
the initial nucleus, p' '^ Pf, this condition is usually fulfilled in the experiments, 
and 

2. the wave function Xp'(k') of the ejected nucleon is negligible outside of an interval 
of momentum Ak' around the central value p'. 

In such cases, we can write 

X^(k')ak'|$i^^)^0. (12) 

This condition is equivalent to assume that the wave function Xp' is orthogonal to the 
components of the initial state. Non-orthogonality effects have been found to be small 
in the region of low missing momentum. However, the assumptions done in this approxi- 
mation are not valid for high missing momentum, for which the region neglected by the 
approximation (|I2|) is explored. 



We are interested here in the low missing momentum region so we can write the current 
matrix element as 

(<l>i^-i),Xp'|>(q)|$S^^) = /rf^A:x;,(k + q)J^(q + k,k)(<|.(;^-i)|ak|$i-'^ (13) 

In this equation we identify the single-particle overlap function between the states $o 
and ^l^~^\ defined, in momentum space, as the matrix element |jl| 

^a(k) = (<l>(^-i)|ak|«l>i^V. (14) 

Using this definition, we can write the many-body matrix element of the current ([TsD as 
a matrix element between single-particle states, namely between the overlap function and 
the distorted wave of the final proton 

(<^L^-'^Xp'l>(q)l$i^^) = {Xr>'\n^)\^a). (15) 

This is the matrix element that we compute in the present work in order to obtain the 
(e, e'p) response functions. The information about short-range correlations is contained 
inside the overlap functions "^a, which are obtained from a correlated OBDM by the 
asymptotic method explained in the next section. The matrix elements (|T5|) are computed 
by performing a multipole expansion of the current operators in terms of Coulomb, electric 
and magnetic operators, and of the outgoing wave function Xp' in partial waves. The 
corresponding response functions (H) are expanded in Legendre functions of cos 6*' (the 
angle between p' and q), and their expressions are given in appendix A. 

The physical interpretation of the overlap function is clear by writing it in the form 

^„(k) = JK<PM (16) 



where (pa is the overlap function normahzed to the unity, and it is usually identified with 
the effective single-particle wave function of the "shell" occupied by the ejected nucleon. 
The spectroscopic factor Sa = {^a\^a) is the occupancy probability of the shell. 

The plane wave impulse approximation (PWIA) will be useful for the physical in- 
terpretation of the correlation effects shown below in terms of overlap functions. In 
PWIA since FSI is neglected, the wave function of the ejected proton is a plane wave, 
Xp'(k) = (5(k — p'), and hence eq. (0) becomes 

(*L^-'\xp'l>(q)l*^^^) = ^'^(P',P)^.(P), (17) 

where we have introduced the missing momentum, p = p' — q, identified with the mo- 
mentum of the proton before the interaction. As a consequence of the above factorization 
property, the hadronic tensor (|^) is proportional to the momentum distribution |\I/q,(p)P 
of the overlap function 

where i(;'^'^(p', p) is the hadronic tensor for a single nucleon with initial momentum p and 
final momentum p': 

i/;^'^(p',p) = J'^(p',p)V'^(p',p). (19) 

In the same way, the response functions are also proportional to the momentum distribu- 
tion: 

W^ = ^^(p',p)|§„(p)f (20) 

where i/7^(p', p) are the response functions for electron scattering by a single proton with 
momentum p. 

When the FSI is turned on, the above factorization is not true anymore but the general 
behavior of the response functions is preserved. The mean effect of the FSI is a reduction 
of the responses due to the absorptive part of the optical potential. We will see below 
that the effects of the short-range correlations are decoupled from the FSI effects. 



3 Correlated model of OBDM and overlap functions 

3.1 Overlap functions 

The basic quantities of interest for our calculations are the overlap functions between 
nuclear states with A and A — 1 nucleons, eq. (p^. We work in coordinate space, where 
the overlap function is 

^a(x) = ($L^-^)|a(x)|$(^)). (21) 



Here x = (r, Sz, t^) is a generalized coordinate including spin and isospin, r is the relative 
coordinate respect to the center of mass of the residual nucleus $ 



(A-l) 
/ 



and a(x) is the 



destruction operator of a nucleon at the point x. We assume that the initial nucleus is in 
the ground state $o with energy Eq , while the residual nucleus remains in an arbitrary 
state $^^"^\ with energy E^~'^\ 



Using the Schrodinger equation verified by the initial and final nuclear states, a system 
of integro-differential equations for the overlap functions can be written EH]. However 



in the procedure explained below to compute these functions for the low-energy levels of 
the residual nucleus we only make use of its asymptotic behavior, which is based on the 
following equation verified by the overlap functions at large distances, r ^ oo 






(r) + (A - l)t;(r)vl/,(r) = [eJ^^ - Ei^-'^]^^iv 



(22) 



where v{r) is the A^A^ potential and r = ri — r2 is the relative coordinate. This equation 
means that the overlap function behaves asymptotically as a single-particle interacting 
with the A — 1 nucleons of the residual system as if they were located at the same position, 
namely at the center of mass of the residual nucleus. Of course this is only valid for so 
large distances that it is not possible to take notice of the small separation distances of 
nucleons within the nucleus. 

In the cases when the initial nucleus has spin zero and the parity of nuclear states 
is a good quantum number, it is possible to separate the overlap function in radial and 
spin-angular parts [^ 

'i^a{r) = (pnij{r)yijmiO,<P). (23) 

with j = Ja the spin of the final nuclear state ^^~^\ and I = j ± 1/2, depending on the 
parity of this state. The radial overlap function (pniji^) (where the quantum number n is 
reminiscent of the later identification with shell-model states) verifies for large distances a 
radial equation coming from eq. (p2D, with an asymptotic eigenvalue given as the difference 
between the energies of the initial and residual nuclei Eq — E^~^\ which is a negative 
number for every value of the excitation energy of the residual nucleus. Therefore the 
overlap function behaves as a bound state and it has the typical exponential decay 



-kr 



c- 



oo 



(24) 



with 



\ 



2/i 



E'c^ 



(A-l) 



-K 



(A) 



(25) 



The exponential decay is modified by a logarithmic phase in the case of proton emission, 
where the Coulomb potential plays a role. However this fact does not modify the following 
results. 



The relation between the overlap functions and the OBDM follows from the definition 
of the density matrix of the initial nucleus: 

p(ri,r2) = E(<V^(ri,s)a(r2,.)|«l>(^)) (26) 

s 

Inserting a complete set of states 1$^"^"^^) of the residual nucleus between the two Fermi 
operators we obtain an expansion of the OBDM in terms of overlap functions 

p(ri,r2)=E^l(ri)^.(r2) (27) 

a 

In this work we consider the OBDM expanded in multipole densities with angular mo- 
mentum / 

p(ri,r2) = ^5]pKn,r2)P/(cos0i2) (28) 

where 612 is the angle between ri and r2. Inserting the expression ( p3D into (^) and 
performing the sums over third components we find the corresponding expansion of the 
OBDM multipoles in terms of radial overlap functions 

Pi{ri,r2) = X](2j + l)0n/j(ri)0„«j(r2) =J2^nijiri)^nijir2) (29) 

nj nj 

where we have defined ipnij{f) as the radial part of the overlap function normalized with 
a factor y/2j + 1. 

iJni,{r) = ^2j + 10„,,(r). (30) 



3.2 Asymptotic methods for computing the overlap functions 

We recall here the method presented in [O to compute the overlap functions by means 



of the exact OBDM for the ground state of the A-particle system. 

We consider a fixed value of /, and denote the corresponding radial overlap functions 
with angular momentum I simply as '(/'^(r), a = 0, 1, 2, . . ., corresponding to bound states 
of the residual nucleus ordered by increasing energy Eq < El < ■ ■ -, and with 
asymptotic behavior 

p rCotCL 

IpaM ~ Ca 5 a ^ 00 (31) 

a 
with K = [2p(^i-^~^) - 4^^)]^/V/i. 

The asymptotic behavior of the OBDM is, from eq. 



p f^aCL 

pi{r,a) r^^^a{r)Ca , a -^ 00. (32) 



Now, due to the ordering ko < ki < ■ ■ ■, the above sum is dominated by the first term, 
with the slowest exponential decay, for long distances a such that a{ki — /cq) ^ 1 



g-fcoa 



pi{r, a) ~ ^poir)Co , a ^ oo (33) 

a 

Combining this equation with the asymptotic behavior of the diagonal part, which allows 
us to determine the constant Co, 



-2koa 



pi{a, a) r^ \Co\ — 5—, a ^00, (34) 



a 



we can compute the overlap function ipoif) and the corresponding separation energy. 
In order to obtain the second overlap function, we apply the procedure to the density 
obtained subtracting the contribution of the first overlap 



p—kia 



piir.a) -i)Q{r)^IJo{a) r-u i)i{r)Ci , a ^ 00. (35) 

a 

In principle all the overlap functions corresponding to bound states of the residual nucleus 
may be obtained by repeating these steps. In the next section we check the validity of 
this procedure, which in the present paper we call "exponential decay method" . 



There is an alternative, equivalent way [15] of obtaining the overlap functions from 



the asymptotic behavior of the OBDM, without using explicitly the exponential decay 
property. We will illustrate it using the fact that the diagonal part of the OBDM has also 
an exponential decay given by eq. (^4|) , from where we can write 



p—koa 



p;(a,a) ~ |Co| , a ^ 00, (36) 



which is precisely the behavior of the first overlap function, eq. (|3lD . Using this equation 
in the asymptotic form of the OBDM, we obtain 

V'o(r) = hm , (37) 

Note that in the case Co < we obtain a minus sign which is just a global phase that can 
be inserted in the overlap function. The application of this expression, which we call "^^^ 
method" , has clear advantages over the exponential decay one when it is used to compute 
the overlap functions from a model OBDM without the correct asymptotic behavior (for 
instance, constructed with harmonic oscillator single-particle wave functions). 



3.3 Model of correlated OBDM 

In this work we compute the overlap functions of closed shell nuclei by applying the last 
method explained to a correlated OBDM. We use the model of ref. |2^ which includes 



10 



short-range correlations up to first order in a cluster expansion of the OBDM. The density 
and momentum distributions were compared with the FHNC calculation of ref. [^, with 
a good agreement between both models. 

We begin with the OBDM of the initial nucleus |$q ), written as 

A r 

p{Xi,X2) = (^) (^) dx2, ...,dxA^i'^^*{xi, ...,XA)^i^\xi, ...,xa) (38) 

(<l>o |$o ) -^ 

For the applications to (e, e'p) reactions, we will only need the proton density, this is 
obtained inserting in the previous equation the projection operator Q{1) = |(1 + ^2(1)). 

Short-range correlations are introduced within the model by the Jastrow ansatz for 
the nuclear wave function 

$(^Hl,2,...,A) = F(l,2,...,yl)4f(l,2,...,yl). (39) 

Here ^^i is a Slater determinant and F is a correlation function containing two-body 
central correlations 

F{i,...,A)= n fin,), (40) 

i>j=i 
where r^ = |rj — rj\, and the function f{rij) has a Gaussian functional dependence 

/(r) = l-Aexp(-5r2) (41) 

We use the parameters A = 0.7 and B = 2.2 fm~^ which were fixed in [Q by minimizing 
the nuclear binding energies for the Afnan and Tang S3 interaction. In ref. |2^ up to six 
spin-isospin correlation channels were included. However, the numerical effort grows in 
the present case, since we have computed the OBDM up to distances so large as 100 fm, 
essential to separate the first overlap function in some cases. This fact compelled us to 
reduce the number of correlation channels and to use a Gaussian dependence, in order to 
perform analytically the multipole expansion of the correlation function /(r). 

The OBDM is calculated by a cluster expansion, writing the correlation function as 

Fil,...,A) = l + l[hin,) (42) 

i>j 

and performing an expansion up to second order in h. The resulting OBDM for protons 
can then be written as 

p?(ri,r;) = f/,{r„r[)+A{r,,r[)-B{r,,r[)-C{r,,r[) + D{r,,r[) 
= Po(ri, r'l) + po(ri, r'l) J rfVaif (ri, r'^, r2)po(r2, r2) 
-/dV.p„(r,.r.)i/(r,.r;.r.)p„(r..r;) 
- /*V, / dSp„(r„ r,)p„(.,. r;)p„(r„ r3)ff (r„ r„ r,) 
+ / <i'r. / <i'r3P„(r„ r,)p„(r„ r^r,. r;)i/(r.. r.. r,) (43) 
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Figure 2: Diagrams considered in the cluster expansion of the OBDM. The dashed lines indicate 
the dynamical correlations f{rij) and the solid lines the uncorrelated density. 
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where the new correlation function H can be expressed in terms of the correlation function 

f{r) 

H{r,,rj,rk) = Q{l)[f{nk)f{rjk) - 1]. (44) 

This function contains the correlations between the two pairs of particles {ik) and {jk). 
The operator Q{1) guarantees that the particle 1 is a proton. 

The functions A, B, C, D corresponding to the corrections to the uncorrelated OBDM, 
p°(ri,r'^), are represented diagrammatically in Fig. 2. Therein, the open circles represent 
the two coordinates ri and r'^^, while the solid dots refer to coordinates of inner nucleons 
which are correlated to the rest. A continuous line represents a non correlated OBDM 
Po(ri, r2), while the dashed lines joins the different particles involved in the new correlation 
function H. Thus in the diagrams A and B, the particles 1 and 1' are simultaneously 
correlated to a third particle 2. On the other hand, in diagrams C and D there are two 
inner particles 2 and 3 which are correlated between them. 

Using the above expression a multipole expansion is performed to obtain the radial 
densities, pi{r, r'), needed to compute the overlap functions for different angular momenta. 
We refer to P2| for details on this expansion. Note that, on the contrary to ref. W^ 



in this model we do not separate the multipoles pij (r, r') of the density in spin-orbit 
partners j = I ± 1/2 explicitly. However these two contributions are included in pi. Since 
our correlated OBDM is based in a single-particle basis obtained with a Woods-Saxon 
potential including spin-orbit interaction, the two overlap functions ipnij corresponding to 
an occupied shell have different energies from the beginning and, in principle, they can 
be separated in the asymptotic region. 



4 Test of the asymptotic methods in the shell model 

Before computing the overlap functions using the correlated model of sect. 3.3, it is 
convenient to perform a test of the asymptotic methods using a nuclear model where the 
exact solution is known a priori. In this way we will be able to determine (i) which of the 
algorithms introduced below is more adequate to extract the overlap functions, and (ii) 
the asymptotic distance needed to separate the several overlap functions. 

We perform this analysis in the extreme shell model (SM), where the overlap functions 
are just the single-particle wave functions of the occupied shells. Apart from its simplicity, 
another reason to choose the SM is that corresponds to the zero-order of the correlated 
model. Under the assumption that the Jastrow correlations, as a first-order correction to 
the SM, do not drastically change the asymptotic behavior of the OBDM, we expect that 
the convergence conditions found in the SM will keep valid in the correlated model. 

The single-particle wave functions in the SM are obtained by solving the Schrodinger 
equation for a nucleon in a Woods-Saxon potential 



Vws{r) = --— 773RT -l-cr 



"0 



l + exp(-^— ^) \m.„c ) r dr 
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Vu 



l + exp(^) 



+ Vc{r), (45) 



Nucleus 




Vo [Mev 


Vis [Mev 


R [fill 


ao[fm 


ais [fm 


12(3 


P 

N 


62.00 
60.00 


3.20 
3.15 


2.74 
2.74 


0.57 
0.57 


0.57 
0.57 


160 


P 

N 


52.50 
52.50 


7.00 
6.54 


3.02 
3.02 


0.53 
0.53 


0.53 
0.53 


40Ca 


P 

N 


57.50 
55.00 


11.11 
8.50 


4.10 
4.10 


0.53 
0.53 


0.53 
0.53 


48Ca 


P 

N 


59.50 
50.00 


8.37 
7.54 


4.36 
4.36 


0.53 
0.53 


0.53 
0.53 



Table 1: Parameters of the Woods-Saxon potential. 

where Vc{r) is the Coulomb potential for protons, of an uniform charge distribution with 
radius R. The parameters of the potential are given in table 1 for the closed shell nuclei 
considered in this work. 

We have solved numerically the radial equation up to distances of r = 100 fm, in order 
to compute the corresponding OBDM in a wide asymptotic region, where we will be able 
to check the convergence of the methods. First, we compared our wave functions with 
the ones obtained integrating the equation up to 11 fm, as is done traditionally, obtaining 
essentially the same answer in both cases up to the region close to r ~ 11 fm. The energy 
eigenvalues of the proton shells are shown in the third column of table 2. For each value 
of /, we construct the shell model OBDM as a sum of the corresponding single-particle 
radial densities of the occupied states with angular momentum equal to /. 



4.1 Exponential decay methods. 

We first focus on the asymptotic decay method, in which the overlap functions are obtained 
by fitting the exponential decay of the OBDM, eq. (B^). This can be done in several ways: 



(I) Logarithm fit. Taking the logarithm in both sides of eqs. (|33| , |3^ ) we have asymp- 
totically. 



\n{a\pi{r,a)\} 
ln{a^Pi(a, a)} 



\n {Co\tlJo{r)\} -koa 
2kQa. 



\nC', 



(46) 
(47) 



We first compute Cq by fitting a straight line to ln{a^pi(a, a)} (we assume Cq > 
since this is just a global phase). Then |'?/'o(^) I is computed by fitting another straight 
line to ln{a|p/(r, a)\} and dividing by Cq. Finally, the sign of ipo{^) is obtained from 
eq. (P^D as the one of pi{r, a). 
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Nucleus 


nlj 


\E^s [Mev 


12(^ 


lSl/2 


32.27 




lP3/2 


15.49 


160 


lSl/2 


27.36 




lP3/2 


13.92 




lPl/2 


9.29 


40Ca 


lSl/2 


35.54 




2Si/2 


9.80 




lP3/2 


26.12 




lPl/2 


22.63 




14/2 


15.83 




14/2 


8.37 


48Ca 


lSl/2 


39.25 




2Si/2 


14.73 




lP3/2 


30.15 




lPl/2 


28.00 




14/2 


19.92 




14/2 


15.14 



Table 2: Energies of the proton single particles in the shell model, Ewg. 

(II) Trace minimization. The overlap function is calculated dividing the OBDM by 
an exponential 

Pi{r,a) 



V'o(r) 



Cn 



(4J 



with Co and ko determined from the diagonal part pi{a,a) as in method I. The 
remaining parameter, a, is chosen by imposing that the overlap density p^{r,r') = 
'^o{'^)'^o{r') be as close as possible to the OBDM, pi{r, r'), for every value of r and r' 
contained in an asymptotic interval [ai, a„]. This condition is achieved by minimizing 
the trace functional 

2 



F{a)^Tr[pi-pl 



dr 



dr' 



ai 



pi{r,r') - 



pi{r,a)pi{r',a) 

/^2cxp(— 2fcoa) 

^0 7? 



(49) 



(III) Trace minimization with three parameters. This third method is a variation 
of fit II, that considers Cq and ko as additional parameters in the trace functional. 
So the three parameters are now fixed by computing the absolute minimum of this 
functional. 



The fit procedure II is similar to the one applied in ref. |jl9[. In the following we 
compare the results provided by these methods for different choices of the asymptotic 
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interval [ai,a,i] where the fit is performed. We shall study the case of ^°Ca in the shell 
model, where the OBDM contains / = 0,1,2 multipoles. 

In Fig. 3 we show results for / = 0, corresponding to the overlap functions of the shells 
2si/2 (first overlap function) and lsi/2 (second overlap). Panel (a) shows an example of 
what is obtained using a OBDM computed up to a^j = 11 fm. In addition we have used 
ai = 2.3 fm, which is the point where the density po{r, r) reaches the 10% of its maximum. 
This value is not large enough to be considered asymptotic and the resulting overlap 
functions are clearly incorrect. In this figure we note a misbehavior of fit I (dashed lines) 
in the region close to the node, where eq. (|33D is not valid, since the OBDM is dominated 
there by the second overlap function. This misbehavior is not found in fits II (dot-dashed 
lines) and III (dotted lines) because in both cases the exponential fit is done globally and 
not point by point. Since the first overlap has not been adequately extracted, we also 
obtain an incorrect result for the second overlap, shown in panel (b). We note in this 
panel that the displayed curves stop around 2-2.5 fm. The reason is that the subtracted 
diagonal part Po{r,r) — 0o(^)^ becomes negative in this region as a consequence of the 
incorrect value of 0o(^)- 

The results for the first overlap function improve when we increase a„ to 15 fm and 
ai to 5.9 fm, corresponding this last distance to the point where the density is 1% of its 
maximum value, as it is shown in panel (c). Even though there is a clear improvement 
with respect to the results of panel (a), there is still a small difference with respect to the 
exact result (solid line), which is larger for the results of fit I. These small differences are 
amplified when the second overlap function, shown in panel (d) , is computed. Nevertheless 
there is also a clear improvement respect to the former results of (b). 

In order to find a reasonable agreement with the exact result we have to use ai = 13.7 
fm, where the density reaches the 10~^% of its maximum value. The corresponding results 
are shown in panel (e), where we have used again a^ = 15 fm. Although not seen in the 
scale of the figure, the results of fits II and III are closer to the exact result than the 
corresponding to fit I. Finally, the second overlap function is shown in panel (f), where 
we still note small differences with the exact result specially for low r. The results for 
the second overlap rely heavily on the adequacy of the computed first overlap function. 
These small differences can be further minimized if a higher value for the asymptotic 
points ai,au is utilized in computing the first overlap. We will see this when we discuss 
the y/p method. 

In the case of the second overlap for / = we use a different interval [a'l, a'^] from 
the one considered in the corresponding first overlap. The upper limit a'^ is chosen as 
the point where the subtracted density pQ{r,r) — ^o(^)^ becomes negative, since this is 
a clear indication that the first overlap is incorrect at this point. The lower limit a'l is 
chosen as the point where the subtracted density is 10% of its maximum value. This 
value is not critical in the cases in which the first overlap function is incorrect, as there 
is no improvement by changing a^. On the other hand, in the cases in which we obtain 
a reasonable result for the first overlap function, the result for the second one is already 
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Figure 3: Overlap functions of '^^Ca for / = in shell model, computed with the exponential 
decay method I (dashed lines), II (dot-dashed lines) and III (dotted lines). With solid lines we 
show the radial functions, corresponding to the exact result. The several panels refer to different 
asymptotic intervals [a/,au] considered in the fit of the first overlap function, namely [2.3 fm, 11 
fm] for (a,b), [5.9 fm, 15 fm] for (c,d), and [13.7 fm, 15 fm] for (e,f). 
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quite good for low values of a^. This is a consequence of the simplicity of the shell model 
where we are working since in this case the subtracted density is factorizable as a product 
of single-particle wave functions of the lsi/2 shell 

Po(r, r') - Mr)Mr') = 2Ru,,,{r)R^s,/,{r'). (50) 

and it is not necessary to separate a third overlap function. Of course in the correlated 
case, where there are extra contributions to the OBDM, one should be careful and choose 
a value of a'l for which there is convergence. 

In Fig. 4 we show the first overlap functions of ^°Ca for / = 1,2, corresponding to 
the shells lpi/2 (a) and ld^/2 (b). The exact results are shown with solid lines, while 
with dashed, dot-dashed and dotted lines we show respectively the results of the fit II 
performed for three different asymptotic intervals with upper limits a„ = 11, 15 and 20 
fm, and lower limits corresponding to the points where the OBDM reaches the 10%, 0.1% 
and 10~^% of its maximum value. In panel (a), corresponding to the lpi/2 shell, the three 
fits give a similar result which is around a factor of two higher than the exact answer. The 
results do not show a noticeable improvement when the asymptotic interval is increased 
from [5.9,15] fm to [13,20] fm. In panel (b), the results of fit II corresponding to the 1^3/2 
shell are again above the exact answer, although we note a convergence trend in going 
from the interval [6.6, 15] (dot-dashed) to [17.2,20] fm (dotted). The results obtained with 
fits I and III are not shown in the figure; fit III gives essentially the same result as fit II, 
while fit I is worse than fit II. 

It is clear from these results, that it is not possible to extract the first overlap function 
for the p and d shells, using asymptotic distances up to 20 fm. This is related to the fact 
that the energies of the spin-orbit partners (lpi/2, IP3/2), and (1^3/2, l'^5/2) have close 
values, their difference being around 3.5 MeV and 7.5 MeV respectively (see table 2). 
Hence, at 20 fm the contribution of the second overlap function is still important (see 
panels (e)-(f) in the same figure), specially in the case of the p shell where distances close 
to 100 fm must be used in order to get convergence, as we will discuss when we study 
the y/p method. Since in all these cases it was not possible to extract reasonable values 
for the first overlap function, we do not show the incorrect results for the second overlap 
function. 

Recently Gaidorov et al. [^ have presented results for the overlap functions of the 
p-shell in ^^O using a method similar to fit II. The procedure was applied to the correlated 
OBDM of [^ computed up to 11 fm, which uses the same shell model we are considering 



here, with similar energies and wave functions for the lpi/2 and lp3/2- An averaged value 
for the Ip wave function was extracted in that reference; we will explore this possibility 
in the shell model. In order to compute an averaged overlap function we first assume that 
the multipoles of the OBDM (|29|) can be approximated by 

pi{ri,r2) ^1^2(2/ + l)0nz(ri)0„z(r2) = 5Z^r.z(ri)i/'„z(r2) (51) 
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Figure 4: First overlap functions of ^'^Ca for / = 1 and / = 2 in shell model, computed with 
the exponential decay method IL With solid lines we show the radial functions of the Ipi /2 and 
1^3/2 shells, corresponding to the exact results. The asymptotic intervals [a;,au] considered in 
the fit for / = 1 are [4.0 fm, 11 fm] (dashed), [5.9 fm, 15 fm] (dot-dashed) and [13.0 fm, 20 fm] 
(dotted); for / = 2 they are [4.5 fm, 11 fm] (dashed), [6.6 fm, 15 fm] (dot-dashed) and [17.2 fm, 
20 fm] (dotted). In panels (c) and (d) the overlap functions have been "averaged" dividing by 
[2(2^ + 1)]^'^ instead of [2j + 1]^'^, and also the wave functions of the 1^3/2 and 1^5/2 shells are 
shown with short-dashed lines. In panels (e) and (f) the asymptotic behavior of the p and d 
wave functions has been displayed. 
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where (f)ni{i^) is a mean value of the two spin-orbit partners (pnijif^), with j = / ± 1/2. The 
overlap function ^niii") is now normalized with a factor [2(2/ + 1)]^^^, and we assume an 
asymptotic exponential behavior ipni{i^) ~ C exp{—kr)/r. Hence we have for the OBDM 

g-fcr' 

pi{r,r') ~ ipni{r)C — —, r' -^ oo. (52) 

Now we can proceed as before, by fitting an exponential decay to the OBDM, and com- 
puting the averaged overlap function as (f)ni{i") = i'ni{i")/[2{2l + l)]"*^^^. 

Results for the averaged overlap functions of the Ip and Id shells of '^'^Ca obtained 
by this procedure are shown in panels (c)-(d) of Fig. 4. The meaning of the lines and 
the asymptotic intervals used in the fits are the same as in panels (a)-(b), but here we 
also include, for comparison, the radial wave functions of the lp3/2 and 1^5/2 shells with 
short-dashed lines. Note that the overlap functions displayed in panels (c)-(d) are related 
to the ones of (a)-(b) just by a global factor [(2j + l)/2(2/ + 1)]^/^^ 

We begin discussing the results for the Ip shell shown in panel (c). The fits shown 
with dashed and short-dashed lines are on average in the intermediate region between the 
lpi/2 and lp3/2 curves. However for higher values of the asymptotic interval (dotted lines) 
the fit begin to move out of this region, — it is now similar to the 1^3/2 wave function — 
and the possibility of obtaining an average value breaks down. The case of the Id shell 
(d) is more intriguing. While dashed and short-dashed lines are between the ^3/2 and ^5/2 
curves, the dotted lines are well below them. The conclusion extracted from these results 
is that only for low values of the asymptotic interval [ai, a„] used in the fit a mean value 
of the overlap function is provided. However the results are unstable, since they change 
when another interval in the fit is used, and depend on the particular /-shell. 

This behavior can be easily explained in the shell model. For example in the case of 
the p shell, the exact OBDM is computed as 

pi(r,r') = 2i?ip^/2(^)^iPi/2(^') +4i?ip3/2(r)i?ip3/2(^)- (53) 

As seen in panel (e) of Fig. 4, the wave functions of the two partners lpi/2 and 1^3/2 are 
quite similar for low r. We note that they are also similar in the asymptotic region up to 
~ 11 fm, since their energies are very close. So we obtain the following approximation for 
the OBDM up to ~ 11 fm 

pi(r,a) ~ 6i?ip(r)i?ip(a) (54) 

and for a large, but not exceeding ~ 11 fm the two wave functions have a similar expo- 
nential behavior so we can write 

Pi(r, a) ~ v6-Rip(r)C , a ^ 11/m (55) 

and then it is possible to extract Rip{r) between the two partners, as it is shown in panel 
(c) with dashed lines, where a < 11 fm. For larger values of a the two wave functions 
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Figure 5: Averaged overlap function of the p shell of ^^O in shell model, computed using fit II 
in the asymptotic intervals [3.5,11] fm (dotted lines), [9.5,15] fm (dot-dashed lines), and [12.2, 
20] fm (short-dashed lines). The exact functions of the lpi/2 and 1^3/2 shells are shown with 
solid and long-dashed lines respectively. 

begin to separate due to the different exponential decay. For the interval [6.6, 15] fm (dot- 
dashed hues in panel (c)) the two wave functions are still quite close and eq. ( ^5] ) remain 
approximately valid. However, for the dotted lines the asymptotic interval is [13,20] fm, 
where the two wave functions are clearly different and eq. (^) is not valid. In this case 
the exact asymptotic behavior 



-k\a 



Pi[r,a) 



4R 



lP3/2 



2R 



lPl/5 



Rip{r) 



[r)C, 



ACie-''^'' + 2006-^°" 



ir)Co- 



-koa 



(56) 



should be used. The exponential fit of this equation may be performed, but the result 
will not be the searched quantity -Rip(r) and will depend on the interval [a/, a^] used and 
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on the fitting method. In addition, in the very far limit where the second overlap can 
be neglected, the fit procedure will converge to the exact first overlap function ipQ{r) = 
^2j + lRip_^^^{r) = \/2Rip_^^^{r). Since the pro- mediated overlap function is computed 

dividing to [2(2/ + 1)]^''^ = V^, by this procedure the "averaged" overlap function will 



converge to the wrong result 0(r) = Jl/SRip.^^^. 

Finally, in the / = 2 case, the energy difference between the 1(^3/2 and 1^5/2 shells is 
bigger than in the case of the p-shell. This makes that the two wave functions separate 
at shorter distances ~ 8 fm, as seen in Fig. 4, panel (f). This implies that the results 
are less satisfactory than for the p-shell for high asymptotic interval (see dot-dashed lines 
in panel (d)). These are clear indications of the impossibility of extracting an averaged 
overlap function using this method. 

The interesting example of the p-shell in ^^O is shown in Fig. 5. Therein we show 
with solid and dashed lines the corresponding wave functions of the lpi/2 and 1^3/2 shells, 
which are very close below 5 fm, where they begin to separate. In the upper panel we also 
show with dotted lines the averaged overlap function obtained with fit II in the asymptotic 
region [3.5,11] fm, which should correspond to the mentioned calculation by Gaidorov et 
al., in ref. |^. The dotted line is between the two p-shell wave functions, but if the same 
fit is applied to higher asymptotic intervals we obtain the dot-dashed and short-dashed 
curves of the figure, that are below the exact result. 



This example clarifies why the results of ref. ||2T| are not far wrong in this particular 
case of ^^O, since they used an OBDM computed up to 11 fm. However the way in which 
the average is done is not under control. So in the correlated case one should be careful 
in the interpretation of the results, because the exact ones are not known a priori, and 
the effects due to correlations cannot be unambiguously separated from the fit procedure. 



4.2 ^ method. 

Next we examine the alternative method based in eq. (|37|) , where the overlap function is 
directly computed as the quotient 

p{r,a) 
Mn = I ^ , (57) 

for a value of a large enough to reach convergence, which can be easily checked by com- 
puting for several values of a. This method has several advantages over the exponential 
fits previously analyzed. First it has not adjustable parameters and no numerical mini- 
mization has to be carried out. Second, the OBDM has not to be calculated in an interval 
[ai,a„], but only in a few asymptotic points a. This is preferable in the correlated case, 
where the computation of OBDM becomes longer. Moreover, when the density is factor- 
izable p{r,r') = ip{r)ip{r'), then eq. (0) always provides the exact overlap function for 
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Figure 6: Overlap functions of ^''Ca in shell model computed with the y^ method. All the 
functions are normalized with a factor \/(2j + l)/47r. In each panel we use the three indicated 
values of the asymptotic point a, that correspond in ascending order to the dashed, dot-dashed 
and doted lines respectively. The exact result is shown with solid lines. The second overlap 
functions (panels on the right) have been computed using the first overlap function from the 
corresponding left panel and a second asymptotic point a' = a — 1 fm. 
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any value of a. This makes this method exact in the shell model for all the shells in ^^C 
and 1^0. 

In Fig. 6 we show results for the non trivial case of ^''Ca in which two overlap functions 
are present for each /. In this case the asymptotic expression of the OBDM can be written 



Piyr.a) 



g-feoa 



V^o(r)Co + ^o(r)Cie('=«-'=^H. (5. 



and a should be chosen large enough in order to neglect the second overlap function. An 

appropriate value of a can be obtained by imposing exp[(fco — ki)a\ = 10~^. This makes 

the second overlap function contribution to be of the order 0.1%. The estimated value of 

a so obtained is 

6 9 

« ^ I— IT- (59) 

However, in practice the convergence will be reached at a different value, depending on 
the relative value of the constants Cq and Ci. In the shell model we will determine the 
adequate value of a by comparing with the exact result. 

The case of the 2si/2 shell of ^°Ca is shown in panel (a) of Fig. 6. Therein we represent 



the first overlap function normalized with a factor J{2j + l)/47r, and computed from 
eq. (|57|) for three different values of a. Namely we show the results for a = 4 f m with 
dashed lines, a = 7 f m with dot-dashed lines, and a = 12 fm with dotted lines. Using this 
last value we already reproduce, within the scale of the figure, the exact overlap function 
shown with solid lines. This value is in agreement with a ~ 11.5 fm provided by eq. 



For each one of the curves presented in panel (a) we compute the second overlap 
function using an asymptotic point a' that must be less than a. This is clear if we 
remember that the method matches the asymptotic behaviors of the OBDM and of the 
first overlap function for distances r > a. Then when we build the subtracted density the 
asymptotic contribution will vanish in this region. This is equivalent to say that the point 
a effectively acts as the infinite point a ~ cxo, so in this numerical method it has no sense 
to compute the second overlap function for a' > a. For this reason, the second overlap 
functions for / = displayed in panel (b) of Fig. 6, have been computed for asymptotic 
point a' = a — 1 fm. Again for the biggest value shown, the extracted overlap function 
almost coincides with the exact result shown with solid lines. 

Results for the remaining shells / = 1,2 are shown in panels (c)-(f). As before, we 
show in each panel three curves corresponding to three ascending values of a indicated in 
the figure, with dashed, dot-dashed and dotted lines respectively. In the case of the second 
overlap function (panels on the right) we use the corresponding first overlap function of 
the left panel and asymptotic point a' = a — 1 fm. The values of a for which convergence 
of the first overlap function is obtained are shown in table 3. In the case of the p-shell 
(panels (c) and (d)), we find convergence for a ~ 100 fm. Using the single-particle energies 
of table 2 and eq.(^), we obtain a ~ 86 fm for the p- shell and a ~ 60 for the d-shell. 
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Nucleus 2gi/2 lpi/2 lc?3/2 

1*^0 80 

4°Ca 12 100 61 

^^Ca 12 > 100 98 



Table 3: Values of the asymptotic points for which convergence of the first overlap function is 
reached with the y^ method. 

We have performed the same study of the proton overlap functions using the ^/p 
method for other closed-shell nuclei. The results of the convergence values for ^^O and 
^^Ca are summarized in table 3. In the case of ^^O we only show the non-trivial case of 
the p-shell. The corresponding overlap functions can be separated by using the OBDM 
computed up to a ~ 80 fm. Finally, In the case of '^^Ca, the separation energy within the 
p- and (i-shells is smaller than in ^°Ca (see table 2). As a consequence, the convergence 
values are larger than in the former case. In particular, for a = 100 fm the lpi/2 and lp3/2 
are not completely separated and the convergence value is not shown in the table. The 
estimated value of convergence for this shell is a = 138 fm. 

As a summary of this section, with our present study using the shell model, we have 
shown the reliability of the asymptotic methods to compute the overlap functions of 
nuclei from the knowledge of the /-multipoles of the OBDM. Our results have shown the 
necessity of studying the convergence of the results in each case and that in many of them 
one should calculate the OBDM up to such huge distances as 100 fm in order to separate 
the overlap functions. In relation to the several methods studied, all of them provide 
the correct result if the asymptotic interval [a;,a„] is within the region of convergence. 
However, due to its simplicity, the y/p method is preferable in the general case in which 
the OBDM is the solution of a correlated many-body problem, and this is the method we 
will use in the next section to compute the correlated overlap functions. 



5 Results for (e, e'p) observables and overlap functions 
in the correlated case 

In this section we present results for overlap functions, spectroscopic factors, (e, e'p) re- 
sponse functions, and cross sections, using the correlated model introduced in section 3. 
Thus we go beyond the single-particle model and will be able to identify the effects of 
short-range correlations on these quantities and observables. 
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5.1 Quasi- hole overlap functions 

We start with the correlated OBDM of closed-shell nuclei, defined by eq. (^3|), and com- 
pute the multipoles, p/(r, r'), as shown in ref. ||2^. The zero-order density p^{r,r') in 
eq. (^) corresponds to the non correlated shell model of sect. 4. The correlated overlap 
functions for quasi-hole states can be obtained by using the y/p method discussed in the 
last section. We apply eq. {^7\) to the correlated OBDM for asymptotic points, a, large 
enough to reach convergence. In the present case the exact result are not known a priori, 
but we are guided by the former study performed in the shell model. It is expected that 
the values of the convergence asymptotic points in the correlated case do not change too 
much respect to the shell model ones. This can be understood in our model by studying 
the correlated OBDM, eq. (|43|) . Using the fact that the correlation function /(r) —>■ 1 for 
r — > CX3, we have for the function H defined in (H4|) 



if (ri, r'l, r,) ~ Q(l) [/(ns) - 1], r[ ^ oo. (60) 

On the other hand, the non-correlated density is dominated asymptotically by the first 
overlap function of the shell model 

Po(ri,r'i) ~ 0o(ri)0o(ri) r[ ^ oo. (61) 

From these equations it is straightforward to obtain the following asymptotic expression 
for the correlated density (^3|) for r[ —>■ oo 

p{{r,,r[) r^ Kir^)Mr'l) r[ ^ oo (62) 

where the function K{ri) is defined as 

Kin) = 0o(ri) + 0o(ri)y'rfV2g(l)[/(ri2)-l]po(r2,r2) 

- J d^r,po{r,,r,)Q{l)[f{ru) - l]0o(r2) 

-/.3,^/.3,3^o(r„r2)0o(r2)po(r3,r3)ii(r2,r2,r3) 

+ J d^r2 J dV3po(ri, r2)po(r2, r3)0o(r3)ii(r2, r2, r3) (63) 

We can also determine the asymptotic behavior of this function K{ri) by using the par- 
ticular Gaussian form (^) of the correlation function, so for ri — » oo the second and third 
terms in (^) can be neglected with respect to the other terms due to its Gaussian decay. 
Then we can write 

K{ri) ~ ?7(/)o(ri), ri ^ oo (64) 

where the constant rj is defined as 



V = 

+ 



1- f d^r2 /"rf3r3|0o(r2)|Vo(r3,r3)ii(r2,r2,r3) 
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Using now the y^ method, the first correlated overlap function reads 

K(r)0o(rO K(r) 

V^o(r) = hm , = ^^. (66) 

'' /ir(r')0o(r') v^ 



r'— >oo 



Since the function K{r) ~ r70o(r), for r ^ oo, we see that this correlated overlap function 
behaves asymptotically as the non-correlated one multiplied by the constant ^ 

^o(r) ~ 0700 (r), r ^ oo (67) 

Hence in the present model the short-range correlations do not modify the energy of the 
first overlap function respect to the shell model, since the asymptotic behavior of the 
OBDM is determined by the exponential decay of the single-particle wave function 0o(r). 
Note that for shorter distances the above proportionality ( |67D does not hold because, in 
that case, the function K{r) includes other terms depending on the non correlated density 
and on the correlation function /(r), as can be seen in eq. (|63[) . 



The same procedure can be applied to each one of the multipoles pi{r, r') to show that 
the energy of the first overlap function for each value of / does not change respect to the 
uncorrelated case. The same conclusion was also obtained in ref. [|15| in a model similar 
to ours by starting with the OBDM pij{r,r'). However in our model it is not possible to 
prove easily a similar result for the energy of the second overlap function. 

These arguments indicate that the asymptotic points needed to compute the overlap 
functions are similar to the ones found in the shell model. In any case, in our calcula- 
tions we have checked numerically the convergence of the results for the different overlap 
functions. For this reason, the correlated OBDM has been computed for values of the 
asymptotic point as high as 100 fm in order to separate the first overlap function in the 
cases in which the energies of two overlap functions are close in the shell model. The 
values of the asymptotic point a where convergence is reached are given in the third col- 
umn of table 4. We first note that the convergence values of the asymptotic point a are 
similar to the ones obtained in the shell model (compare with table 3). Thereby, in order 
to extract the lpi/2 overlap function, we need to go up to ~ 86 fm for ^^O and up to 
~ 100 fm for '^"Ca, while in the case of the 1^3/2 overlap function for ^°Ca, convergence 
is found for a ~ 64 fm. This indicates that in fact the separation energies of correlated 
overlap functions are close to the ones of the shell model. 

The results for the correlated overlap functions of the nuclei ^^C, ^^O and "^"Ca are 
shown in figs. 7-9 respectively. In each one of these figures we show with solid lines the 
radial density r'^ip{rY of the correlated overlap function and with dashed lines the non- 
correlated result corresponding to the shell model. We do not show the overlap functions 
of the also studied nucleus ^^Ca. Having stopped our calculation at 100 fm, it was not 
possible to obtain convergence in this nucleus for the overlap function of the p-shell, which 
requires higher values of a. 
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Figure 7: Solid lines: radial density of overlap functions for ^^C computed from the correlated 
OBDM. With dashed lines we show the non-correlated results in the shell model. 

In figure 7 we show the first overlap functions of ^^C for / = 0, 1. In both cases 
convergence is reached for relatively small values of a = 7 and 8 fm, respectively, due to 
the fact that in the shell model there are no second overlap functions contributing to the 
OBDM. Short range correlations introduce extra contributions in the OBDM. However, 
these extra contributions decay much faster than the exponential one, and very large 
distances are not needed to extract the first overlap function. As we can see in figure 
7, the inclusion of short-range correlations produces in both cases a reduction of the 
maximum of the overlap function in coordinate space, while there is an increase for high 
r, which is better seen in the case of the p-shell, since it lies at higher distances (panel 
(b)). In this last case we observe in addition that the overlap function undergoes a small 
shift to the right due to the correlations. 

Similar effects are observed in figs. 8 for ^^O and fig. 9 for ^"^Ca. In all cases there is 
a reduction of the overlap function at intermediate distances (in most of them coinciding 
with the maximum of the radial density) and an increase for more large distances. As 
in ^^C, we also observe a shift to the right of the overlap functions corresponding to the 
outer shells. This is the case of the shells lpi/2 and lp^/2 in ^^O and 2si/2j lc^3/2 and ld^/2 



m 



40 



Ca. 



This shift effect over the outer shells can be understood in terms of the repulsive 
nature of the NN interaction for short distances, implicit in the correlation function /(r), 
and the well known healing property of the wave function for the two-nucleon system. 
The correlation function produces a wound in the A^A^ wave function \I', and what we are 
seeing in the overlap function is the average effect of healing due to the interaction of the 
outer nucleons with the nucleons in the core. 

However the inner shells do not show this effect because the short-range repulsion 
due to the core partially cancels the one produced by the external shells. The net effect 
depends on the particular nucleus and on the shell involved. For instance, in the case of 
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Figure 8: The same as fig. 7 for tlie nucleus O. 

the lsi/2 shell in ^°Ca, shown in panel (b) of fig. 9, the correlations produce a shift of the 
overlap function to the left, i.e., into the nucleus, since the short-range repulsion by the 
outer shells tends to compress the Is wave. The same compression effect is observed in 
the inner lobe of the 2si/2 (see panel (a) of fig. 9). 

In the cases of the intermediate shells lpi/2 and lp3/2 for '^''Ca the joint effect of 
repulsion by the inner and outer shells produces a shift to the left in the low r region and 
a shift to the right for large r. Hence the net effect of correlations over these shells is a 
widening of the overlap function, as seen in panels (c) and (d) of fig. 9. 

Once the overlap functions have been extracted, we can compute the corresponding 
separation energy by a fit to a function (Je'^^'"^^^^^ /r for large distances. We have 
performed this fit in the interval between 11 and 28 fm for the correlated and uncorrelated 
overlap functions obtaining essentially the same energies. The inclusion of the logarithm 
Coulomb phase is important in this fit for protons, since it can modify the extracted 
energies in more than 2 MeV in the case of ^°Ca. The relative difference between correlated 
and uncorrelated energies is shown in the fourth column of table 4. In all cases the 
differences are less than 0.5 % even in the case of the ^^Ca shells where convergence was 
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Figure 9: The same as fig. 7 for tlie nucleus ^''Ca. 
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Nucleus nlj a [fm] \Ec — Ewsl/E^js [%] Sth 5', 



cxp 



i^C lsi/2 7 4 X 10-5 0.985 0.59 

lp3/2 8 < 10-^ 0.986 0.56 

1^0 lsi/2 7 2 X 10-5 0.985 

lp3/2 7 X 10-4 0.986 0.59 

lpi/2 86 4 X 10-5 0.986 0.57 

4°Ca lsi/2 6 X 10-2 0.988 0.75 

2si/2 14 2 X 10-3 0.992 0.64 

lp3/2 3 X 10-2 0.985 0.72 

lpi/2 100 3 X 10-3 0.985 0.72 

1^5/2 2 X 10-3 0.985 0.74 

1^3/2 64 3 X 10-5 0.985 0.74 

4^Ca Ts^2 6 X 10-4 0.986 
2si/2 16 1 X 10-^ 0.991 
lp3/2 2 X 10-^ 0.946 
lpi/2 > 100 3 X 10-^ 1.058 
1^5/2 1 X 10-3 0.983 
1^3/2 100 6 X 10-4 0.983 



Table 4: For each one of the quasi-hole states in the closed shell nuclei studied we show: the 
asymptotic distance a for which convergence of the correlated overlap function is reached (third 
column), the relative difference between correlated and uncorrelated separation energies (fourth 
column), and the computed spectroscopic factor (fifth column). For comparison we show also 
the experimental value of the spectroscopic factors extracted from (e, e'p) experiments. 

still not found for the overlap function. These numerical results confirm that short-range 
correlations do not change the mean field values of the separation energies for quasi-hole 

states. 

Results for the spectroscopic factors resulting from our model are shown in the fifth 
column of table 4. These have been computed as the norm of the corresponding correlated 
overlap function 

S=m)- (68) 

As seen in table 4, all of the computed spectroscopic factors are slightly less than one, 
being in most of the cases around S ~ 0.985. The only exception found in table 4 is the 
value Sp^ = 1.058 for ^^Ca, which is not a definitive number since the asymptotic point 
a = 100 fm used is not large enough to reach convergence in this case. 

Our results indicate that short-range correlations of Jastrow-type reduce the shell 
model occupation probability no more than 2%. This reduction is not enough to explain 
the experimental values extracted from (e, e'p) analysis shown in column 5 of table 4. This 
is in agreement with other studies which report values similar to ours for the spectroscopic 
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factors. Van Neck et al. found in []I5| that central correlations generate a reduction 



around 1-2% of the occupancy probability in ^^O. More recently Fabrocini and Co' |^| 



have computed overlap functions within the FHNC/SOC theory and report values around 
0.97-0.99 for the spectroscopic factors with Jastrow correlations. Spin-isospin and tensor 
correlations (not included in our calculation) produce additional reduction of these values 
to 5 ~ 0.86-0.9 for the valence shells. The discrepancy with experimental values could be 
further reduced by the inclusion of long-range correlations [jl^. Center of mass corrections 
however produce an enhancement of ~ 7% of the p-shell in ^^O |^. Further investigation 
including all of these effects in a consistent way is needed in order to clarify the situation. 



5.2 Quasi-particles and the continuum 

Up to now we have restricted our study to overlap functions corresponding to quasi-hole 
states. Our correlated OBDM model allows us to compute also the multipoles pz(r, r') for 
high values of /, which are expected to contain contributions coming from quasi-particle 
states, i.e., states non occupied in the shell model but which are partially populated in 
the ground state due to nuclear correlations. We have investigated if the overlap functions 
for quasi-particles can be extracted from our OBDM using the asymptotic method. 

This study is motivated by a recent calculation done in [|l^] where results for the 



quasi-particle overlap function for the Irf-shell in ^^O and for the 1/ shell in ^''Ca are 
presented. These authors begin with an OBDM including Jastrow correlations and apply 
an asymptotic procedure similar to fit II in order to extract the overlap functions. For 
instance they report a value of 5* = 0.01 for the spectroscopic factor of the 1/ shell in 
^°Ca. Neither the asymptotic interval used for the fit nor the convergence distance are 
indicated in [^ . Apparently they should not have used very high asymptotic values since 
they use a harmonic oscillator basis that fall off rapidly at large distances. 

However these results were criticized in ref. [|^, where it was shown that it is not 
possible, starting from a CBF-type wave function, to generate bound-state overlap func- 
tions with quantum numbers that are unoccupied in the Slater determinant. The reason 
is that the overlap functions decay exponentially with the same decay constant as the 
hole single-particle orbital 

In fact, when we apply the ^p method to compute a quasi-particle overlap function 
from our radial OBDM for high values of / we do not obtain convergence within the 100 
fm range and instead the results decrease rapidly, being negligible for high r. We must 
conclude that in our model it is not possible to obtain such quasi-particle states. 

This result can be understood by examining the following asymptotic expression of 
the OBDM for unoccupied multipole / 

Piifur'^) r^ Cf{ri,r[), n^r'^^co (69) 



32 







Figure 10: The solid lines are the correlated OBDM, p2{r,r') for / = 2 in ^^C, as a function of 
r for different values of r' = 6, 7, 8, 9, 10 y 11 fm. The dashed lines are a fit to an asymptotic 
approximation f{r,r') defined in (^0|). 



where C is a constant and the function f{r,r') is given by 



f{rur[) 



exp {—2kri) exp {—2kr[) exp ( — yi]_liL\ 



rfr'i'^lri + r'lY 



(70) 



Here k is the wave number of the (occupied) valence shell and B is the parameter of the 
correlation function. This expression is proved in Appendix B for the simplest case of the 
multipole / = 2 for the nucleus ^^C 

In figure 10 we show with solid lines the computed radial density p2{r,r') for ^^C as 
a function of r for several fixed values of r'. In addition we show with dashed lines the 
function /(r, r') multiphed by a convenient constant C fitted to the density. We see that 
in fact the above asymptotic expression is approximately verified by the computed density 
for high values of r and r'. 

If we now try to compute a quasi-particle overlap function using the y/p method we 

obtain 

pi{r,a) 2a , , _Rr^_„^2, 



^/(r) 



pi[a,a) 



r'^{r + a) 



;Po[r)e 



-B{r-aY/2 



(71) 



this expression goes to zero for a — ;> oo, since the diagonal part of pi has an exponential 
decay, while the non-diagonal part has an additional Gaussian behavior corresponding to 
the correlation function. This explains why in our results the extracted overlap function 
for quasi-particles are zero. Thus in our model of correlated OBDM it is not possible 
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to obtain the overlap functions for quasi-particles since the corresponding information of 
single-particle states above the Fermi level (configuration mixing) has not been included 
into the model. Moreover, when one subtract from the correlated OBDM the contribution 
of the quasi-hole states the remaining density contains only the contributions coming 
from the continuum states of the residual nucleus. These contributions are implicit in the 
expansion (^) in terms of overlap functions and they can be expressed as an integral over 
the energy. 

p(r, r') - Pguas^-hole{T, t') = J rfE^s(r)t*s(r') (72) 

With the asymptotic method studied here it is not possible to extract these overlap 
functions of the continuum. For this a practical inversion method of the integral ( [72D in 
the asymptotic region is needed. The knowledge of these overlap functions would be of 
interest, for instance, to compute the (e, e'p) cross section for high missing energy. 

5.3 Exclusive response functions and cross sections 

In figs. 11-13 we show the five exclusive responses for proton knockout from the valence 
shells of the nuclei ^^C, ^^O and ^°Ca, as a function of the missing momentum. In all 
the cases the kinematics correspond to a fixed value of the momentum transfer q = 460 
MeV/c and u fixed around the quasi-elastic peak. In each panel of figs. 11-13 we show 
four curves corresponding to different models for the initial and/or final wave functions 
that enter in the current matrix element (|^). We show results for PWIA and DWIA with 
and without short-range correlations in the initial state overlap functions. The DWIA 
results have been obtained using for the FSI the optical potential of ref . |^ . Specifically, 



the solid lines have been computed with the DWIA model and using correlated overlap 
functions, while the dashed lines do not include correlations. Thus comparison between 
solid and dashed lines shows the effect of short-range correlations in the responses. Results 
in PWIA with and without correlations are shown with short-dashed and dotted lines, 
respectively. 

In all the cases we note in the region p < 200 MeV/c an increase of the L, T and TL 
responses due to correlations, which is around 5% near the maximum. This increase is 
quite independent on the FSI since it is also present in PWIA. The reason for this fact 
is that correlations between the ejected proton and the residual nucleus have not been 
included. The increase seen in the responses can be easily understood in PWIA, where 
we are basically seen the momentum distribution of the shell, as a consequence of the 
hardening effect seen in the overlap functions in figs. 7-9 for the valence shells. Since in 
momentum space the low p region is sensitive to the high r region, the increase of the 
overlap function for high r translates into an increase of the Fourier transform for low 
p (~ 100 MeV/c) where the maximum of the momentum distribution is located. Since 
the correlations in the ground state are in some degree decoupled from the FSI, the same 
effect is propagated to the case of the DWIA. 
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Figure 11: Response functions for proton knock-out from the P3/2 shell in C, for q = 460 
MeV and lo at the quasi-elastic peak. Results are shown in DWIA with (solid lines) and without 
(dashed lines) correlations in the hole overlap function, and in PWIA with (short-dashed lines) 
and without (dotted lines) correlations. 
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Figure 12: Response functions for proton knock-out from the valence shells, 1^3/2 and lpi/2, in 
^^O, for q = 460 MeV and u at the quasi-elastic peak. The meaning of the lines is the same as 
in fig. 11. 
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In the case of the TT response, we also find an increase of its absolute value due to 
correlations in DWIA, which is not seen in PWIA because the leading-order magnetization 
current do not contribute to this response |^ , and the resulting factorized single- nucleon 
TT response is of order (p/M)^ in a non-relativistic expansion (this is the reason why this 
response function is so small). This kinematical dependence comes exclusively from the 
convection current, producing a hardening of the maximum of the momentum distribution 
toward higher p-values ~ 150 MeV, where the correlated and uncorrelated results are 
closer. In DWIA, the FSI breaks the factorization property and the magnetization current 
gives a contribution which wherefore is much larger than the PWIA result. 

Regarding the fifth response function TL', which only can be measured using polarized 
electrons, it is exactly zero in absence of FSI. In DWIA however it produces a contribution 
to the total cross section and the correlations in the ground state produce an increase 
which is in general of the same order as in the unpolarized responses. This increase is 
even more large (~ 15%) in the case of the 1^3/2 shell in ^°Ca (see fig. 13). 

State-independent short-range correlations produce an increase of the (e, e'p) cross 
section, since the later is a linear combination of the several response functions appearing 
in eq. (1). An example is shown in fig. 14, where results of DWIA calculations for 
the reaction ^^0(e,e'p) are displayed together with the experimental data of ref. [^ 



Here the kinematics correspond to fixed momentum transfer q = 570 MeV/c and energy 
transfer cj = 170 MeV, at the quasi-elastic peak. The energy of the electron beam is 580 
MeV and the proton is emitted in the scattering plane. In fig. 14 we show with solid lines 
our DWIA results using the correlated overlap functions for the lpi/2 and 1^3/2 shells, 
while with dashed lines we show the uncorrelated results. We note an enhancement of the 
cross section due to correlations which is of the same order of magnitude as was found for 
the response functions, and which clearly increases the disagreement between theory and 
experiment. In the same figure we show with short dashed and dotted lines the computed 
cross sections multiplied by the factors 0.6 (lpi/2) and 0.5 (IP3/2) in the correlated case 
and 0.64 (lpi/2) and 0.53 (1^3/2) in the uncorrelated one. 

Hence the scaling factor needed to reproduce the experimental cross section is smaller 
for correlated than for uncorrelated overlap functions even though the computed spectro- 
scopic factors for these shells are smaller than one, S = 0.985 (see table 4). This fact does 
not necessarily imply a decrease of the experimental spectroscopic factors since these are 
obtained by a simultaneous fit of the parameters of the single particle potential also. In 
other words, experimentally one searches for the best phenomenological overlap function 
that when included into an uncorrelated DWIA code reproduces the experimental data. 
Our model has not adjustable parameters since the correlations are already included and 
so are the spectroscopic factors. Our results are showing that short-range correlations of 
the central-type in the ground state do not produce an improvement of the (e, e'p) data 
description. Tensor correlations and long-range correlations are obvious candidates for a 
reconciliation between theory and experiment. 

To end the discussion we give in fig. 15 another application of our DWIA model in 
photo-nuclear reactions. Therein we show the computed (7,p) cross section from the 
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Figure 13: Response functions for proton knoclc-out from the valence shells, 14/2 and 14/2) 
in ^'^Ca, for q = 460 MeV and io at the quasi-elastic peak. The meaning of the lines is the same 
as in fig. 11. 
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Figure 14: Computed (e, e'p) cross section for tlie valence sliells of ^^O. The solid lines include 
correlated overlap functions while the dashed lines do not. These calculations have been scaled 
to give the dashed and dotted lines in order to reproduce the experimental data from ref. 



lpi/2 shell in ^^O for two beam energies of E^ = 60 and 72 MeV together with the 
experimental data of ref. [^ . Again we show with solid and dashed lines the correlated 
and uncorrelated results respectively. No scaling factors are included. The impact of 
central correlations in this case is completely different from the (e, e'p) reaction. In this 
case they produce a large reduction of the cross section. The difference between the two 
reactions lies in the different kinematical region which is being probed by photons. In the 
case of (7,p) the energy-momentum transfer verifies u = q and we are far from the quasi- 
elastic peak region. As a consequence the missing momentum is well above 200 MeV/c. 
Thus apart from the small values of q, the photon is exploring here the high momentum 
tail of the overlap function where correlations produce a reduction of the momentum 
distribution. Such reduction can also be appreciated in figs. 11-13, where the correlated 
transverse response (the one contributing to photo-reactions) is below the uncorrelated 
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Figure 15: Computed {'j,p) cross section for the lpi/2 shell of O. The solid lines include 
correlated overlap functions while the dashed lines do not. Experimental data are from ref. 



one for p > 200 MeV/c. One should also be aware of the difficulties presents in the DWIA 
description of the (7,p) reaction for so low energies and high missing momenta, where in 



particular the orthogonality approximation (10) is no longer true and other effects [^ of 
the same order as the ones arising from correlations could appreciably change the results 
of fie. 15. 



6 Summary and Conclusions 

In this work the effects of short-range correlations on (e, e'p) observables and overlap 
functions have been studied. The starting point for the present calculation has been the 
OBDM, including short-range correlations of central Jastrow type, and which has been 
computed by a cluster expansion to leading order in the correlation function. Correlated 
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overlap functions corresponding to quasi-hole states in closed-shell nuclei have been ex- 
tracted from the asymptotic OBDM multipoles, pi{r,r'), computed at large asymptotic 
distances r' < 100 fm. The reliability of the extraction method has been tested in the 
shell model, where a detailed study of the different fit procedures and of the convergence 
distances has been performed. These distances have been found to be very large in the 
cases in which two single-particle states of close energies are present. In those cases the 
asymptotic distances considered in our calculation of the OBDM have been enough to 
separate the corresponding overlap functions of the nuclei ^^C, ^^O and "^^Ca. As also has 
been found in other works |]22| , [33| , ^], short-range correlations produce small effects on 
the density distribution and likewise on the OBDM. Our results given in figs. 7-9 show 
that these effects are more noticeable in the overlap functions since they are determined 
by the asymptotic behavior of the OBDM. Short-range correlations produce a redistribu- 
tion of the single-particle densities in coordinate space. Their values are reduced at the 
maximum and increased for large distances. In the case of the valence shells we find a 
hardening of the distribution related to the collective effect of the NN repulsion at short 
distances. 

The values of the computed spectroscopic factors in the present work are around 
0.985, in accord with the findings obtained with other techniques. It is known that tensor 
and long-range correlations can reduce these values but, at present, no model is able to 
reproduce the experimental values extracted from (e, e'p) data. 

The computed overlap functions have been included in a DWIA model of the (e, e'p) 
reaction, and exclusive response functions and cross sections have been computed for 
quasi-elastic kinematics. Although the computed spectroscopic factors are less that one, 
we have found an increase (~ 5%) of the response functions and, accordingly, of the 
cross section in the region of the maximum of the missing momentum distribution for 
knockout from the valence shells. This reduction is independent on the FSI and is a 
consequence of the increase of the single-particle densities for large distances. Thus the 
inclusion of central correlations is worsening the description of the experimental data in 
our model. This again proves that central correlations alone are not enough to describe 
this reaction successfully. Apart from spin-isospin and tensor correlations, not included 
here for computational reasons, we would like to remark the necessity of a model including 
in addition long-range correlations in a consistent way. 

The later correlations are related to the presence of admixture of Tcmlti-hu configu- 



rations into the valence wave functions of the residual nucleus [|^ . For instance in the 
case of the ^°Ca nucleus, the residual states correspond to the nucleus ^^K. It is known 
that the transverse form factors of the measured elastic and inelastic transitions in ^^K 
show significant departures from the single-particle picture and that a modification of the 
extreme shell model wave functions through the effect of core polarization is needed to 
describe the electron scattering cross section |^. When computing the (e, e'p) reaction 
in ^°Ca one uses a wave function which reproduces the elastic electron scattering data or 
equivalently, the charge density, corresponding to the initial state (the ground state of 
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^°Ca), ignoring the necessity of a proper description of the residual states also, usually 
treated as single holes in the core. 

Concerning the (7,p) reaction, central correlations play here a more important role, 
producing a large reduction of the cross section, since this reaction is sensitive to the high 
momentum components of the wave function where correlation effects are maximized. 

In this paper we have demonstrated with a realistic model that the asymptotic method 
to compute the OBDM is feasible and that convergence of the results can be obtained: 
therefore it rises as an alternative, reliable starting point to be applied to other kind of 
correlated densities. 
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A Multipole analysis of exclusive response functions 

The general multipole analysis of the (e, e'p) responses including polarization degrees of 
freedom of the target nucleus and electron was presented in refs. p3| , p9| . The formalism 
can also be applied to the present case of J = nuclei, where some simplification of 
the multipoles given in ref. |2^ can be done. The following equations we write in this 



appendix have been obtained, after some work, from the corresponding equations of ref. 
[ P3[] for the particular case JT = {J' being in [^ the angular momentum corresponding 
to a multipole expansion in terms of spherical harmonics Yjm{(^*,4'*) of the nuclear 
polarization angles). 

We expand the nuclear electromagnetic current as a sum of Coulomb (for the time 
component), and electric and magnetic (for the transverse three- vector current) multipole 
operators of rank J. The final hadronic state is also expanded in partial waves of the 
ejected proton, as a combination of hadronic states with total angular momentum J, 
denoted as \(r) = \{lj)Ja] J), that represents a nucleon in the continuum with asymptotic 
angular momenta Ij coupled with the residual nuclear state I'l'^"^"^-') = \Ja)- The exclusive 
response functions can be written in the form 

W^ = ^EW^l(cos^')W^l (73) 

W^ = ^Y.[L]Pli^ose')Wl (74) 
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(75) 
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with [L] = v2L+T. Note that the whole dependence on the emitted proton angle 6' is 
given through the Legendre functions P^^{cos9') and that the present response functions 
are divided by the factor K = Mp'/{27ih)^ respect to those of ref. [|^. The reduced 
response functions W^ , defined as the coefficients in the expansions ([7^-[77|) are given by 
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Note that the TL and TL' reduced response functions of ref. p3[ include an extra factor 
-\/2 due to the different definition of the vtl and vtl' factors. The coupling coefficient 
^cr'a{L) includes the internal sums over third components and it is defined as 
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We also use the parity functions 
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In order to define the functions R^, ^ and I^, ^ in eqs. (ff^)- 
electric and magnetic multipole matrix elements 

a = {a\\Mj{q)\\0) 
E^ = {o\\ff{q)m 
M. = {a\\tfr\q)m. 



(83) 

(84) 

(85) 
(86) 

we introduce the Coulomb, 

(87) 
(88) 
(89) 
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where Mj{q), ff{q) and Tj'^'^^q) are the usual Coulomb, electric and magnetic multipole 
operators. The functions R^, ^ and I^ ^ in eqs. (|78D -(|82D are then defined by the following 
quadratic forms constructed with these multipoles 

(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 
(97) 
(98) 

The L, T, TL and TT responses include only the real parts of the quadratic combinations 
of the various multipole matrix elements, while the fifth response function TL' is a linear 
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combination of the imaginary parts (^,0). Therefore the TV response is zero in PWIA, 
where the matrix elements (p7|- |89|) are real numbers. In presence of an interaction, how- 
ever, the matrix elements (|87|-p9D are in general complex numbers, due to the asymptotic 
complex phase e"^'J introduced by the nuclear interaction in the wave function, and as a 
consequence, the fifth response function is different from zero in DWIA. The sum over the 
quantum numbers a = {l,j, J), a' = {l',j', J') and L in eqs. (|73D-(^2D is only restricted 
by angular momentum conservation. In practical calculations we fix the number of mul- 
tipoles involved in the sums by comparing our results with the ones corresponding to the 
factorized PWIA in the impulse approximation where, as known, the nuclear responses 
can be computed exactly. 

The outgoing proton wave function corresponds to a solution of the Schrodinger equa- 
tion for positive energies using a complex optical potential fitted to elastic proton-nucleus 
scattering data. The partial wave Ij with energy E > and wave number k = \plME is 
determined by the asymptotic condition 



RiAk, r) ~ J^^e-*("'+'''^) sin (kr - j] log2fcr - /- + a^ + b^S (99) 

^ irh k \ 2 -^ J 

where Sij is the complex phase-shift and ai is the Coulomb phase-shift. In the limit 
in which the imaginary part of the optical potential is zero, the phase-shift 6ij a real 
number, the continuum radial wave functions are normalized with a Dirac delta function 
containing the energies (see Ref. ||23|). The imaginary (absorptive) part of the optical 



potential modifies the normalization of the continuum sates since the fiux of the outgoing 
particles in the elastic channel is reduced. 
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B Asymptotic OBDM for unoccupied states 

We consider as an example the simplest case of the multipole / = 2 for the nucleus ^^C 
More details are given in |3^. We use the expression P3| ) for the correlated contributions 
to the OBDM. First we exclude the contribution of diagrams C and D of fig. 2, since 
the dependence of these terms on the density coordinates r, r' is done across the non 
correlated density po(r; r') which only contains the multipoles / = 0, 1 in the case of ^^C 
In other words, the external points r and r' in diagrams C and D are connected with the 
others with density lines only, which cannot modify its multi-polarity / = 0, 1. 

In the case of diagram A of fig. 2, it can be written in a multipole expansion as [^] 

nihji Ih ^ ~^ \ / 

X / rfr2r^po(r2)4(ri,r2)4K,r2) (100) 



where the function fi2{ri, r2) is the multipole of the correlation function for angular mo- 
mentum I2. The sum over ni,/i,ji corresponds to the occupied states lsi/2 e lp3/2 in 
^^C The multipole / = 2 of the OBDM is obtained as the coefficient of the Legendre 
polynomial P2(cos^ii') in the above equation. Since we are interested in the asymptotic 
behavior for ri,r[ ^ cxd we only consider the contribution coming from the h = 1 term, 
i.e., (riiliji) = lp3/2- Hence the 3-j coefficient gives a non zero result for /2 = 1 only. The 
corresponding multipole ^2 = 1 for the correlation function is proportional to the integral 



/rfcos^i2 cos^i2/(ri2) = -A f dcosOuCoseue-''''^ ~ -— g-^^^^-''^)'. (lOi; 



Hence the asymptotic behavior of the / = 2 multipole for ri, r^ -^ 00 is 



00 _ p-B(ri-r2)2 g-B(ri-r2)2 



PA(ri,r^)i=2~i?ip3/2(ri)i?ip3 K) / dr2rlpo{r2)- -, . (102) 

' ' Jo r\r2 T\r2 

Changing to the variable r'^ = T2 — Vm-, where Vm = (^1 + ^i)/2 is the mid point between 
ri and r[, we arrive to 

PA{n,r[)i=2 ~ Rip,,M)RiP.;M f '' ? Poirm) H dr',e-'^^'-\ (103) 

TiTi J —00 

Finally, introducing the asymptotic behavior of the radial functions 

Rip,,M) ~ ^ (104) 
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we obtain the following asymptotic expression 

exp (— 2A;ri) exp (—2kr[) exp ( — ]^ '"^ 



2 

rfr'i^lri + r'l " 



pA{ri,r[)i=2 n,^. , ,^2 ^1. K ^ oo. (106) 



A similar expression can be obtained for diagram B of fig. 2. 
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